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Abstract 

This article discusses the interstellar extinction curve in the visible and the value of the ratio of absolute to selective 
extinction R v = A V /E(B-V). It is concluded that the visible extinction curve is likely to be linear in the visible and that 
indirect estimates of Ry from tentative determinations of Ay or from infrared and UV observations are questionable. 
There is currently no evidence of any variation of Ry with direction. If Ry is close to 3, as it has been inferred from 
mid-infrared data, starlight in the visible is extinguished by a factor F/Fq = (2.5e~ 2,lm l A ) E( - B ~ v \ But if the visible 
wavelength range alone is considered, 4 appears as its most natural and probable value and F/Fq - e - 2E ( B - v V A . 
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1. Introduction 

Between 1937 and 1941 a series of papers specified the characteristics of interstellar extinction in the visible, 
provided a better unders tanding of the processes at work, and set the general frame-work for all forthcoming studies. 



Greenstein 



( 1937 | j938[ l proved that extinction of a star's light in the visible was not of the Rayleigh type (oc e fl ^ 4 ) as 



previously fhoughf] but was very close to a decaying exponential of l/A (oc e~"t A , that is the visible extinction law was 
linear in l/A). Such laws are common in nature (atmospheric aerosols, Fig.[T] are the best but not the only example), 
and attributed to size distributions of particles, with no specific implication on their chemical composition. These 
distributions are usually associated with large albedos and strong forward scattering phase functions, two properties 
which in the case of interstellar dust were confirmed in 1940 (Hen yey~& Greenstein|1940 1941 1. 



Linear extinction laws are defined by two parameters, their slope, proportional to the reddening E(B- V) = Ag—Ay 
on the line of sight, and the extinction at a specific wavelength A^ . The extinction at wavelength A normalized by 
the reddening, A^/E(B - V), depends on the absolute extinction to color excess ratio (Ry in the Johnson photometric 
system) alone (Sects. [2]and[3]). 

Greenstein found good spatial homogeneity of interstellar extinction. Assuming that the extinction law goes to in 
the infrared (Stebbins, Huffer, & Whitford 1939 ) a value of Ry = 4 is obtained (see also Mendoza ( 1965 )). In contrast 



Greenstein & Henyey ( 1941| > found from a statistical analysis of the mean reddening at two different wavelengths that 



interstellar extinction becomes zero at 3.2 //m, which implies Ry = 3.3 (+0.1). 

Since the 1940s the debate on interstellar extinction has focused on four essential questions. How close to a 
perfect exponential is the extinction of starlight in the visible? Does extinction vary from direction to direction (from 
an interstellar cloud to another)? On how many parameters does it depend? How does it constrain the composition of 
interstellar dust? 

The purpose of this paper is to evaluate which logical conclusions on interstellar extinction in the visible and the 
near-infrared (4300 A - 1 .4 fim) may be reached from an analysis of the abundant and often contradictory literature 
on the subject: what is the shape of the visible extinction curve? Is the key parameter R v a constant or is it highly 
variable, as different studies separately claim, and what is its value(s)? Can extinction in the UV or in the infrared be 
related to the visible extinction and to the value of Ryl 

Sect.[2]gives the relevant definitions that are used in this paper. Sects. [3] to preview the main features of interstellar 
extinction in three wavelength domains which have been used to investigate the value of Ry, the visible, the infrared 



'The first hint for a A a , with a close t o 1, wavelength dependence of interstellar e xtinction came from |Hall|jl937} , in opposition to previous 
suggestions of a Rayleigh l/A 4 extinction jTrumpler 1930, Stebbins & Whitford 1936 . 



Preprint submitted to Elsevier 



March 6, 2012 



1 .000 



0.100 



0.010 - 



0.001 




.6 



1 .{ 



2.0 



2.2 



2.4 



Figure 1: Night occultation of Sirius observed by satellite GOMOS. Bottom spectrum is the spectrum of Sirius observed at a tangent altitude of 
14.8 km, normalized by the spectrum of Sirius outside the atmosphere. Top spectrum is the same after correction for ozone and Rayleigh (nitrogen) 



extinctions (JVo, = 3.7 10 20 c 



Nn 



1.2 10 cm ). Only aerosol extinction remains which behaves as an exponential of 1/A (red trace, 



oc e ' ). Courtesy of Jean-Paul Bertaux and Laurent Blanot (Service d'Aeronomie, Verrieres-le-Buisson, France). 

(near and mid-infrared, 0.8/vm < A < 3/im), and the UV. Sect. [6]questions the variability of Ry- Sect. [7] attempts to 
give the most logical responses to the questions raised in the previous paragraph. 



2. Definitions 

Let F A be the flux of a reddened star, measured on earth, and Fq A that of an unreddened star of same spectral 
type. The quantity which is readily accessible through observation is the excess of extinction of the reddened star at 
wavelength A over the extinction at a given wavelength Aq, 



E(A-Aq) = -2.5 log 



JLl 

Fq,a 

F ->0 

fo..i„ 



(1) 



Extinction in different directions can be normalized to the same amount of reddening between Ao and another given 
wavelength A\(< Aq) by dividing E(A - Aq) by E(A\ - Aq). Normalized extinction curves are defined by 



E(A - Aq) 
E{A\ - Aq) 

2 



E{A X - Aq) 



R 



(2) 



with 



R 



A 1,^0 



(3) 



E(Ai - A ) 

gAt,A a (MX) measures the extinction per unit amount of interstellar matter along the line of sight. 

Normalized extinction curves g AuAo in different directions all pass through the points (l//lo, 0) and (1 /A\, 1) in the 
(l/A,g Al j ) plane, and have the same mean slope, (A Q l - A^ l )~ l . For extinction to be the same in different directions 
normalized extinction curves need to be identical. In this case knowledge of E(Ai - Aq) fixes the amount of extinction 
and the quantity of interstellar dust on the line of sight at all wavelengths. If g AuAo varies with direction extinction 
depends on E(B-V) and one or more additional parameters. 

Two systems of normalization, (A3, A2) and (A\,Ao) can be converted one into the other by 



gx 3 M (1M) = 



gA^Aoii/A)- g Al ,A (l/A 2 ) 
gA u A (Mh) ~ gAuAa^lh) 



(4) 



In Nandy's system (|Nandy||l964t l/Ap = 2.22 /zm _1 and 1 /Ai = 
unity. In the Johnson & Morgan ( 1953 1 photometric system (1/Aq 



1.22 /mi 1 were chosen to give a mean slope of 
= l/A v = 1.82 /mr 1 , l/Ai = l/A B = 2.27 fmT 1 ) 



the slope is 2.22, and normalized extinction curves are defined by 



with 



god /A) 
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(5) 



(6) 



cr B and cry are the extinction cross-sections of interstellar dust in the B and V bands in the direction of the observation. 

Eq.[5]emphasizes the importance of R v . If R v is known, A A becomes measurable at all wavelengths. Eq.[6]shows 
that Ry depends only on the extinction properties of interstellar dust in the B and V bands (along the line of sight). 
Identical normalized extinction curves in different directions certainly imply equality of Ry in these directions. 

The constancy or not of Ry and its exact value(s) determine our ability to evaluate A A in any direction, to correct 
observations of far-away objects for foreground reddening, and to parametrize interstellar dust models. The Schlegel, 
Finkbeiner, & Davis ( 1998]) derivation of A A in different photometric bands from E(B - V) for instance relies on the 



assumption that Ry is constant and equal to 3.1; their estimate of interstellar extinction in any direction is widely used 
today. 

There is no direct and straightforward way to estimate Ry. E(B - V) is deduced from a star spectral type and flux 
ratio between the B and V bands but Ay can not be ascertained unless the star absolute luminosity and distance are 
precisely known. Existing determinations of Ry rely on indirect methods, the extension of go to the infrared (Sect.Q, 
estimates of A v , or on the UV extinction curve (Sect.[6]l. 

3. The visible extinction curve 



With the exception of Divan's work ( |Divan|1954 i spectral studies of interstellar extinction in the visible are rare. 
There is a lack of a large, reliable data-base of stellar spectra free from atmospheric extinction, such as exists in the UV 
(the IUE data-base). Spectra in analytical form did not become common before the 1950-1960s and a precise removal 
of atmospheric extinction (which also has a dependence on l/A, Fig. [TJ is not obvious. Photometry, sometimes from 
only a few data-points, has been the standard means to study extinction in the visible, although it is less accurate than 
spectral analysis can be. 

The upper plot of Figgis the ratio F/Fo, in the [1250 , 8000 ] wavelength range, for the star HD46223^ With an 
excellent precision the visible part of the spectrum follows a decaying exponential of l/A (and diverges from it in the 
UV). 



observations from J.F. Le Borgne library of stellar spectra, http://www.ast.obs-mip.fr/articlel81.html Le Borgne 2003 



3 




Figure 2: Extinction of HD46223, HD62542, HD192273. From top to bottom, ratios of the spectra of HD46223 to the unreddened star of same 
type HD269698 (15 Mon), of HD62542 to HD32630, of HD 192273 to HD31726. The red lines are the exponentials e -2MB-v)/A^ with A(B _ y) the 
B-V difference between the reddened and the reference stars. The extinction of HD192273 is linear throughout the visible-UV wavelength range, 
in the visible to the near-UV for HD62542, and in the visible only for HD46223. The F/Fq ratios have been normalized so that the fits be exactly 
e -2&(B-v)/A (which i mp ii es Ry - ^ Sect.[7j. 
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This exponential decrease of the extinction of starlight in the visible is confirmed by most st udie^j of visible 
normalized extinction curves. [ETivan| ( |1954] l, |Whitford| ( |1958| l, rNandy| ( [T964) l, |Ardeberg & Virdefors| (|T982) and more 
recently Bondar et al. ( 2006[ l found the same gi,,,i (l//i) in all directions (over the visible domain). 

The most extensive data-set of extinction curves in the visible was gathered by Nandy in the 1960-1970s (see 
review, Wickramasinghe ( [ 1 998| l) . The extinction curve weighted over a large number of observations that he derived 
in 1964 ( |Nandy|1964| is reproduced in Fig. [3] It is to a high degree of accuracy linear in the visible from 2.3 fim' 1 
(4350 ) to 0.8 jum -1 (1 .2 /urn) in the near- infrared ( |Wickramasinghe| 1 998| [Whitford| 1 958} . In this wavelength range 
the exponential dependence of F^/Fqj can equally be expressed as 



= C e 



-a/A 



or by 



Aa = 



{A B -Ay)l\ 1 



-(a/A+b) 



+ Av 



(7) 



= 2.22£ (B - V) 



1 fim 



+ 0.46 (R v - 4) 



From Eqs.[7]and|8] 



The extinction per unit reddening is 



2.02E(B - V) 
0.92E(B - V) X (R v 
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(9) 
(10) 



(11) 
(12) 
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4. The infrared extinction curve 

Observed extinction curves are linear only in the visible and part of the near-infrared domain. Fig. [3] shows that 
the curve deviates in the near-UV above l/A ~ 2.3 yurrT 1 (shortly after the fi-band) where there is less extinction 
than expected. The error margin on Fig.[3]increases dramatically because extinction in the UV no longer depends on 
E(B-V) alone (Sect. [5]). 

In the infrared extinction is much less than at shorter wavelengths. Error margins are large and the extinction 
curve is difficult to ascertain. Infrared extinctio n data ar e limited to a few photometric bands, and thermal emission 
from dust grains shows up above a few fim ( | Sherwood ( 1975) and references therein). Below 0.8 finT 1 published 
mid-inf rared extinction laws do not fo llow the visible law ( Stead & Hoare| (p009), figs. 1 and 2 in Whitf ord] ( |1958| l, 



fig. 2 in Ardeberg & Virdefors|(|1982|>). The extinction curve seems to flatten shortward of l/A ~ 0.8/um - I ~(1.2yum) 



and extinction may be larger than predicted by the linear behavior in the visible. 

Mid-infrared extinction was used ( Johnson|1965 1968 Wegner|2003|l to probe variations of Rv (go(l/A — > 0) — > 
Rv, Eq.[5]). These fi ndings are contradicted by several studies (Bondar et al. 2006; Sherwood|1975| Rozis-Saulgeot| 



1956 Schalen|1975 Schultz|1 975 ), sometimes even in the same regions of the sky. There is today no tangible proof 
of spatial variations of the extinction in the infrared. The legitimacy of using the infrared to determine R v is also 
questionable, especially with respect to the spatial constancy of go(A) in the visible wavelength range (Sect. [7] and 
|Mendoza| ( fT965] l). 



two notable exceptions are the photometric observations of Johnson (fig. 12 in |johnson|(l965) ) and of |Schild| j 1977} who finds that 'in no octave 

of the spectrum is the interstellar absorption strictly proportional to reciprocal wavelength it is becoming increasingly clear that variations with 

Galactic longitude are important' (see also Borgman (1954); Johnson & Morgan 1 1955 1; Johnson 119 68)). Johnso n's observations and methods 
(especially his summary in |Johnso"n| j 1 968| ) are criticized in several articles jSherwood|1975||Rozis-Sa ulgeot 1956 Schalen 1975 , S chultz|1975| . 
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5. The UV extinction curve 



The UV extinction curve exhibits the well-known 2200 bump feature (upper plot of Pig . [2]» ; the extinction is less 
than expected from the continuation of the visible extinction curve (Figs. |2j. Early in the 1970s it was recognized 
that normalized go{l/A) curves depend on direction: E(B - V) alone is not enough to determine the whole extinction 
curve. These properties have largely contributed to the development of new interstellar dust models (see |Hoyle &| 
|Wickramasinghe| (| 1 99 1 |l) . 

Cardelli, Clayton, & Mathis (1989) (CCM in the following) proved in 1989 that most normalized curves with a 
bump can be deduced from a function of l/A which depends upon a single parameter, R ccm . In the CCM framework 
the size of the bump depends on E(B - V) only (this is a crude approximation, see Savage (1975)) while the free 
parameter determines the average slope of the extinction curve (with respect to the one in the visible) and its shape 
in the far-UV. Cardelli, Clayton, & Mathis ( 1989 ) assume that R ccm is equal to R v , an affirmation to be discussed in 
Sect. [6] 

Cases also exist where the linear visible extinction extends to the near-UV, eventually to the far-UV (Fig. E}. A 
linear extinction curve over all the visible-UV can be found in the directions of very low column density within the 
Galaxy (Fig. |2j bottom plot, and Zagury (2001)); similar extinction curves are observed in the Magellanic Clouds 
along lines of sight with even higher column density (Zagury 2007). In some directions, HD62542 (Fig. [2} middle 
plot) or HD29647 in our Galaxy, Sk-69228 and Sk-70116 in the LMC ( |Zagury|2007l ), the extinction is linear down 
to the bump region and diverges from linearity in the far-UV only. All these directions have no bump and cannot be 
fitted by the CCM function. The CCM fit can be improved to include linear extinction curves ( Za gury |2007| ). 



6. On the variations of R v 



If, in the visible, normalized extinction curves go(l/ A) do not depend on direction (Sect. [3} there is no reason why 
the relative extinction cross-sections of interstellar grains in the R and B bands (Eq. [6) should vary, and R v should 
remain constant (Sect. [2]). 

Paradoxically the constancy of Ry is disputed not so much from the study of interstellar extinction in the visible 
itself but from indirect and generally more speculative methods, which often involve the extinction in other wavelength 
ranges. They rely either on infrared data (which are not conclusive, Sect.Q; on tentative estimates of Ay; or on the 
CCM fit of UV extinction curves (Sect. [5}. 

Direct determination of A v uses star count methods or open clusters distance estimates. Star counts for the same 
region have resulted to contradictory conclusions regarding the nature and variability of Ry ( |Johnson|196 8 ; Schalen 
|1975 ; Balazs et al. 2004). The values of Ry range from 3 to 6 (Balazs et al. 2004| ). Distance estimates of open clusters 
are no more convincing. In a recent article Turner (201 1 1 shows that the distance to open cluster Shorlin 1 should be 
revised from over 10 kpc to less than 3 kpc and thereby demonstrates the large errors which are associated to distance 
estimates by means of photometry. The same paper also lists some thirty determinations of R v in the Carina region 
which range from 3 to 5.2. Such variations over a small area of the sky are obviously difficult to reconcile with the 
constancy of go found by Nandy, Divan, Bondar, and others (Sect. [3}. 

The same remark applies to the UV. The use (since 1989) of the CCM fit to probe large and local variations of Ry 
is questionable. The free parameter R ccm of an extinction curve is obtained from eq. 1 in CCM 



Ax 

Ay 



a(x) + 



b(x) 

Rrrm 



(14) 



where a(x) and b(x) are functions of x - 
2.27/irrT 1 ) Eq. 14 yields, in any direction 



or 



A~ alone (Cardelli, Clayton, & Mathis 1989 i. For the B band (x B 



A B . . b(x B ) 
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Figure 4: Dotted spectrum is the ratio of HD48099 (reddened star) to the slightly reddened (E(B - V) = 0.07) star HD47839, corrected for 
extinction. If the slight reddening (E(B - V) = 0.07) of HD47839 is considered, the CCM fit is the red dashed spectrum, with free parameter 
Rccm = 1-2. If half only of this reddening is corrected, the CCM fit is the green spectrum, with R ccm = 3.4, close to the value found in Cardelli, 
|Clayton, & Ma this 1 1989 1. A slight error in the adopted value of E(B — V) can thus lead to largely different values of the CCM free parameter. The 
plot also shows that R ccm has no impact on the size of the bump, which, in the CCM framework, is controlled by the reddening E(B - V). 
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If Rccm — Rl 



1 - b(x B ) 
1 - a(x B ) ' 



(15) 



Should R ccm be equal to Ry Eq. 15 shows that it would have to be a constant (from Eq. 15 and eqs. 3a and 3b in 
CCM: Ry = 1.48) independent of direction. This is in contradiction with the fact that the free parameter of the CCM 
parameterization, R ccm , varies according to direction. Ry is therefore unlikely to be the free parameter of the CCM 
function. 

It is further difficult to conceive how Rccm > which in practice fixes the average UV slope of the CCM fit with 
respect to the mean slope of go(lM) in the visible (fig. 4 in Cardelli, Clayton, & Mathis ( 1989] > and Zagury (2007)), 
can influence the visible extinction curve, let alone Ry. It also should be noted that the parameter R ccm is highly 
sensitive to the estimated value of E(B - V) (Fig. ffl, and is not necessarily the same whether the CCM law or its 
improved formulation in Zagury (2007) are employed. Last, R ccm barely agrees with the value of Ry derived by other 
methods, which find, using similar data-sets, either different values in the same region (in the Magellanic Clouds for 
instance (|Zagury|2007]l) or a constant R v (|Bondar et al.|2006[|Schlafly et al.|20T0|. 



7. Discusssion 



The interstellar extinction curve is now accessible through the entire spectrum although studies in each of the three 
wavelength regions (the infrared, the visible, and the UV) are generally done separately and rely on observations with 
different quality, sensitivity, and resolution. In the infrared there is no conclusive evidence that interstellar extinction 
should depend on another parameter in addition to E(B - V) (Sect. |4| and |Sherwood ( 1975[l). 



Extinction in the UV, although the best documented, remains to be understood. The Cardelli, Clayton, & Mathis 



( 1989 1 paper had a strong impact and was a remarkable breakthrough. It proved that interstellar extinction in the UV 
V) and on an additional parameter only. The CCM function, representative of normalized extinction 
is however purely empirical and void of physical significance. Its free parameter is 



depends on E(B - 
curves with a bump at 2200 

unlikely to be R v (Sect.[6]i. The nature of the parameter which governs UV extinction is still to be determined. 

Over sixty years of observations support a linear extinction law in the visible, with gc(l/A) normalized extinction 
curves i ndependent of directi on (Sect. [3]). Linear extinction laws are attributed to grain size distributions (<f>(a)da = 



a da in van de Hulst ( 1981 )). They have no implication for the chemical composition of the grains responsible for 



the extinction: extinction by aerosols in the atmosphere, like interstellar dust, follows a similar law. The extinction at 
a specific wavelength is due to particles of size close to the wavelength. 

The spatial uniformity of go(lM) m the visible strongly favors a constant value of Ry (Sect. [2]). Estimates of Ay 
in different directions, which essentially are of a statistical character (Balazs et al. (2004! and Sect, ffl, suggest that 
Ry should be between 3 and 6. 

Whitford|(l958| l), which interrupts 



The flattening of the extinction curve observed in the near-infrared (fig. 2 in 



the linear decrease of go(l/A) towards the longest wavelengths, was used to derive a value of Ry close to 3. If Ry is 



indeed close to 3, Eq. 13 may be rewritten as 



= 2.5e- 2l,m/A 



(16) 



<W£(fl_y)=i 



It is however arbitrary to determine Ry from infrared data since below 0.8 //nr 1 the interstellar extinction curve seems 
not to follow the same analytical law as in the visible. Ry should depend on the distribution of those interstellar dust 
grains with sizes close to Ag and Ay and be determined from the extension to longer wavelengths of A,\IE(B - V) 
rather than from the observed mid-infrared data. The mid-infrared flattening, should it be confirmed, may as well 
result from a process independent from the extinction at shorter wavelengths, from a change in grain-size distribution 
for the largest interstellar grains for instance. 

A general difficulty in discussing these estimates of Ry is that Ry can be defined in different ways which do not 
necessarily have the same meaning. The exact definition, Ry = Ay/E(B - V), involves only the reddening in the B 
and V bands. Because direct determination of these quantities is particularly difficult alternative definitions rely on 
the limit of go(l/A — > 0), and on the supposition that A A (l/A — > 0) — > which is true only if there is no gray (neutral) 
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extinction. But the go(lM — * 0) limit can also be considered in two ways, either by taking the limit of mid-infrared 
observations, as it is generally done, or by continuing Nandy's linear extinction law (Pig. [3j to the longer wavelengths. 

As stated above the infrared method, which leads to (Rv)m ~ 3 and to Eq.[l6] is not consistent with the original 
definition of R v . If there is a flattening of the extinction in the near-infrared (Rv)ir must be smaller than A V /E(B - V) 
and underestimates Ry. 

If Nandy's curve in the visible alone is considered, provided that the Z?-band lies within the limits of the linear part 
of normalized extinction curves (that is, if Ag is not or little affected by the near-UV departure from linearity, Fig. [3}, 
and assuming no gray extinction, Eqs . [5] and [T2| necessarily imply Ry = 4. This conclusion can be re-formulated in 
the following way. If linear normalized extinction curves in the visible were to be considered alone one would expect, 
from any size distribution of particles, that the coefficient in front of the exponential in Eq.[T6]be less than 1 (instead 
of 2.5 ). In absence of neutral extinction, which is presumed to be the case for interstellar dust ( Dufay|1954 Sherwood 
1975) , this coefficient, C in Eq. |7](or e -0-92(Rv-4) in Eq [T3J, should be 1 . Then 



Fx 



2.2 fim 



E(B 



V) 
Rv 



(17) 

(18) 
(19) 



Rv = 4 is precisely the value |Turner| ( |201 1 J > has obtained from his most recent photometric measurements in the 
Carina region. This value is not "anomalous" as it was often suggested, but rather appears as a logical outcome of the 
linearity of interstellar extinction in the visible. 

This paper was motivated by the large and logically incompatible discrepancies that appear from a survey of the 
existing literature on the determination of Ry- Depending on the study Ry has been found to be either constant and 
independent of the line of sight or highly variable from one direction to another. 

It is my conclusion that the normalized visible extinction law is most likely linear and independent of direction 
so that R v should be the same constant for all directions. Although a value close to 3 is generally preferred, R v = 4 
would more naturally fit with the linear extinction law found in the visible. 
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